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ABSTRACT 

We apply the theory of algebraic polynomials to analytically investigate the transonic proper- 
ties of hydrodynamic accretion onto non-rotating astrophysical black holes. We first construct 
■ the equation describing the space gradient of the dynamical flow velocity of accreting matter. 

Such equation is isomorphic to a first order autonomous dynamical system. Application of the 
Q . fixed point condition enables us to construct an nth degree algebraic equation for the space 

variable along which the flow streamlines are defined to possess certain first integrals of mo- 
tion. The constant coefficients for each term in that equation are functions of certain specified 
initial boundary conditions. Such initial boundary conditions span over a certain domain on 
the real line R - effectively, as individual sub-domain of KxK (spherical flow) and IxlxM 
(accretion disc) for the poly tropic accretion, and of K (spherical flow) and RxR (accretion 
disc) for isothermal accretion. The solution of aforesaid equation would then provide the criti- 
cal (and consequently, the sonic) point r c . The critical points itself are permissible only within 
£■ — . a certain open interval ]r g , L->oo [, where r g is the radius of the event horizon and L^oo is the 

physically acceptable maximally allowed limit on the value of a critical point. 

f*"*- ■ For polynomials of degree n > 4, analytical solutions are not available. We use the 

Sturm's theorem (a corollary of Sylvester's theorem), to construct the Sturm's chain algo- 
rithm, which can be used to calculate the number of real roots (lying within a certain sub- 
domain of K) for a polynomial of any countably finite arbitrarily large integral n, subjected to 
certain sub-domains of constant co-efficients. The problem now reduces to identify the poly- 
nomials in r c with the Sturm's sequence, and to find out the maximum number of physically 
acceptable solution an accretion flow with certain geometric configuration, space-time metric, 
and equation of state can have, and thus to investigate its multi-critical properties completely 
analytically, where the polynomials in r c are of n > 4 (for complete general relativistic ax- 
isymmetric flow, for example, where n — 14), and thus, for which the critical points can not 
be computed analytically. 

Our work, as we believe, has significant importance, because for the first time in the 
literature, we provide a purely analytical method, by applying certain powerful theorem of 
algebraic polynomials in pure mathematics, to check whether certain astrophysical hydrody- 
namic accretion may undergo more than one sonic transitions. Our work can be generalized 
to analytically calculate the maximal number of equilibrium points certain autonomous dy- 
namical system can have in general. 
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1 INTRODUCTION 



Black holes are the vacuum solutions of Einstein's field equations in general relativity. Classically, a black hole is conceived as a singularity 
in space time, censored from the rest of the Universe by a mathematically defined one way surface, the event horizon. In astrophysics, 
black holes are the end points of the gravitational collapse of massive celestial objects. 'Observed' astrophysical black holes may be broadly 
classified into three different categories, the stellar mass (Mbh^ a few Mq), the intermediate mass (significantly more massive than the 
stellar mass black holes but far less massive than the super massive black holes) and super massive (Mbh^W 6 Mq) black holes. All of the 
above mentioned candidates accrete matter from the surroundings, provided that the sources for such infalling material do exist. Depending 
on the intrinsic angular momentum content of accreting material, either spherically symmetric (zero angular momentum flow of matter), 
or axisymmetric (matter flow with non-zero finite angular momentum) flow geometry may be invoked to study an accreting black hole 
system. Since the black holes manifest their presence only gravitationally, and no spectral infor mation can directly be obtained for these 
candidates, one must rely on the accretion processes to understand their observational signature (Pringle 1981; IKato. Fukue & Mineshige 
19981 : iFrahk e7ai]|2002h . 



The local Mach number M of the accreting fluid can be defined as the ratio of the local dynamical flow velocity to the local velocity of 
propagation of the acoustic perturbation embedded inside the accreting matter. The flow will be locally subsonic or supersonic according to 
M(r) < 1 or > 1. The flow is transonic if at any moment it crosses M — 1. At a distance far away from the black hole, accreting material 
almost always remains subsonic (except possibly for the supersonic stellar wind fed accretion) since it possesses negligible dynamical flow 
velocity. On the other hand, the flow velocity will approach the velocity of light c while crossing the event horizon, whereas the maximum 
possible value of sound speed (even for the steepest possible equation of state) would be c/v3, resulting M > 1 close to the event horizon. 
In order to satisfy such inner boundary condition imposed by the event horizon, accretion onto black holes exhibit transonic properties in 
general. 

A sonic/transonic transition in black hole accretion occurs when a subsonic to supersonic or supersonic to subsonic transition takes place 
either continuously (usually from a subsonic to a supersonic transition) or discontinuously (usually from a supersonic to a subsonic transition). 
The particular value of the spatial location where such transition takes place continuously is called a transonic point or a sonic point, 
and where such crossing takes place discontinuously are called shocks or discontinuities. In supersonic black hole accretion, perturbation 
of various kinds may produce shocks, where some dynamical and thermodynamic accretion variables changes discontinuously as such 
shock surfaces are crossed. Certain boundary conditions are to be satisfied across the shock, and according to those conditions, shocks in 
black hole accretion discs are classified into various categories. Such shock waves are quite often generated in supersonic accretion flows 
having small amount of intrinsic angular momentum, resulting the final subsonic state of the flow. This is because the repulsive centrifugal 
potential barrier experienced by such flows is sufficiently strong to brake the infalling motion and a stationary solution could be introduced 
only through a shock. Rotating, transonic astrophysical fluid flows are thus believed to be 'prone' to the shock formation phenomena. 
The study of steady, s tanding, stationary shock w aves produced in black hole accretion and related phenomena thus acquired an important 
status in recent years (lFukuelll98ill98ll2004lll|chakrabartilll989l:lKafatos & Yanalligg^lYang^ & Kafatos|[l99l : ICaditz & Tsurutil998l: 



Fukumura & TsuratJl2004|:lTakahashi et alJl99alDas||2002| : |Das et alfcooil Abraham et all20od:lpas et alfcoollLu et al.lll997l:lLu & 
20041: iNakavama & Fukuel[l9 89: Naga kura & Yamadal2 008 : Nak avamalll996l : lNagakura & YamadJ2009l :lToth. Keppe ns. & Botche\fl99 



Das & Czern 



2009) 



0„ 

•98; 



On the other hand, a physical transonic accretion solutions can mathematically be realized as critical solution on the phase portrait 
(spanned by dynamical flow velocity/Mach number and the radial distance) of the black hole accretion. This is because, from analytical per- 
spect i ve, pr oblems in black hole accretion fall under the general class o f nonli near dynamics dRav & Bhattacharieefcoo'ljAfshordi & Paczvnskil 



20031: liavl 



2003iJlb1:)Rav & BhattachariedF 



2005a b; Chaud hurv etailbood : I5av & Bhattacharieell2006l . l2007al : iBhattachariee & Ray 2007; 



Goswami e t al. 2007; Bhatt achariee et al.ll2009l) . since accretion describes the dynamics of a compressible astrophysical fluid, governed by a 
set of nonlinear differential equations. Such non-linear equations describing the steady, inviscid axisymmetric flow can further be tailored to 
construct a first order autonomous dynamical system. Physical transonic solution in such flows can be represented mathematically as critical 
solutions in the velocity (or Mach number) p hase plane of th e flow - th ey are associated with the critical points (alternatively known as 
the fixed points or the equilibrium points, see Ijordan & Smitr] J 19991) and lchiconj J2006t) for further details about the fixed point analysis 
techniques). To maintain the transonicity such critical points will perforce have to be saddle points, which will enable a solution to pass 
through themselves. 

Hereafter, 'multi-critical' flow refers to the category of the accretion flow configuration which can have multiple critical points accessible 

to the accretion flow. For low angular momentum axisymmetric black hole accretion, it may so happen that the critical features are exhib- 

J 1 J] I 1 1 

ited more than once in the phase portrait of a stationary solution describing s uch flow (Liang & Thompson 1980; Abramowicz & Zurek 1981 

Muc hotrzeb & PaczvhskJl982l;fM uchotrze bl 19831:^^ & Rat. 

Al 



Zjt. L/ LJL< A 11 L 1_ y 11J1V1 ± S UZ. , LI ± S LI ~J f A lyUJ, A V yj I f i-L/ . CI , IjLI A V L> , AVI 11 L^u \— * Z_L 1 11 y lyUV), J 1U1 UlllKJ VV ILt Ut, l_i_llL' 

19891: 1 Abramowicz & Chakrabartilll99d:lKafatos & Yang||l994l : lYang & Kafatoslll995l : ICaditz & Tsurutalll998l : lDasll2002l : lBarai et alJl20o4 



Abraham 



etalJ 



2006 



Das ct : 



al.l l2007l : 



Das & Czernvfl2009r) . and accretion becomes multi-critical. In reality, such weakly rotating sub- 



Keplerian flows are exhibite d in various phy s ical s ituations, such as detached binary systems fed by accretion from OB stellar winds 



i Illarionov & Sunvaevl dl975h : lLiang & Nolarj fl980h ). semi-detached low-mass non-magnetic bi n aries (|Bisikalo et al.l dl998)). and super 



massive black holes fed by accretion from slowly rotating central stellar clusters (IllarionoyJ ( 11988 ): Ho ( 1999) and references therein). Even 
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for a s tandard Keplerian accretion disc, turbulence may produce such low angular momentum flow (see, e.g. Jlgumenshchev & Abramowicz 
(1999), and references therein). 

Two issues are, thus, of utmost importance in studying a shocked accretion flow around astrophysical black holes. Firstly, the inner 
boundary condition near the horizon dictates the presence of a saddle type sonic point through which the flow must pass to become supersonic 
before it finally plunges through the event horizon. Secondly, shocks in supersonic flow forces the flow to become subsonic. Hence, for 
accretion onto black hole, presence of at least two saddle type sonic points are thus a necessary, but not sufficient, condition for the shock 
formation. Since a sonic point is associated with the formation of a critical point in the accretion flow, multi-criticality, thus, plays a crucial 
role in studying the physics of shock formation and related phenomena in connection to the black hole accretion processes. However, 
depending on the symmetry of the flow and t he equation of state o f the accreting material, such critical points may or may not be isomorphic 
to the corresponding sonic points (see, e.g. jDas & Czernv 20091) for further detail about the clear distinction of the sonic points and the 
critical points). Hence the investigation of the multi-transonic shocked accretion flow around astrophysical black holes finally boils down to 
the following set of operations: 

One constructs the corresponding autonomous dynamical systems, then identifies the saddle type critical points of the phase trajectory 
of the flow. Lastly, the global understanding of the flow topologies are performed - which necessitates a complete numerical integration of 
the nonlinear stationary equations describing the velocity phase space behaviour of the flow. 

However, for all the importance of transonic flows, there exists as yet no general mathematical prescription allowing one a direct analyt- 
ical understanding of the nature of the multi-criticality without having to take recourse to the exis ting semi-analytic approach of numerically 
finding out the total number of physically acceptable critical points the accretion flow can have. lDas. Chattopadhvav & Chakrab arti (2001) 
analytically computed the location of the critical points in an axisymmetric accretion, b y being restricted to a spe cific single case where 
only the polytropic accretion was studied under the influence of a specific potential, the Paczvnski & Wiita i 198fJ) pseudo-Schwarzschild 
potential, and did not provide any generalized scheme for predicting the number of critical points an accretion flow can have. In general, 
one needs to introduce a first integral of motion in the form of a polynomial of the critical point, and involving the initial boundary con- 
ditions ._Hieii_one_jriakesatterrrpt to analytically solve that polynomial to find the roots (the location of the critical points). For adiabatic 
flow in IPaczvnski & WiittJ Jl980h potential, such a polynomial is of degree four, and hence can be solved analytically. However, the specific 
choice of the pseudo potential or departure from the pseudo-Newtonian framework may provide a polynomial in the critical point of degree 
higher than four, and a polynomial of degree n > 4 is non-analytically solvable for i ts roots. Such highe r degree polynomials are indeed 
encountered for some specific flow configuration, for pseudo-Newtonian flow under the lNowak & Wagoner! ( 1 99 1 ) potential, or for complete 
general relativistic axisymmetric flow i n the Schwarzschild metric where the polyn omial will be of degree fourteen. Hence analytical solution 
of the polynomial (as has been done in Das. Cha ttopadhva v & Chakrabartil j200lh ) is not a useful method to explore the multi-transonicity 
for realistic flow structures. Instead of explicitly solving the polynomial, one should rather investigate how many physically acceptable roots 
the polynomial can have. 

This is precisely the main achievement of our work presented in this paper. Using the theory of algebraic polynomials, we developed 
a mathematical algorithm capable of finding the number of physically acceptable solution the polynomial can have, for any arbitrary large 
value of n (n is the degree of the polynomial as mentioned earlier). Hence for a specified set of values of the initial boundary conditions, it 
can mathematically be predicted whether the flow will be multi-critical (more than one real physical roots for the polynomial) or not. This 
paper, thus, purports to address that particular issue of investigating the transonicity of a generalized flow structure without encountering the 
usual semi-analytic numerical techniques, and to derive some predictive insights about the qualitative character of the flow, and in relation 
to that, certain physical features of the multi-criticality of the flow will also be addressed. In our work, we would like to develop a complete 
analytical formalism to investigate the critical behaviour of the flow structure for all possible flow configuration around a non rotating black 
hole. Both spherically symmetric as well as axisymmetric flow governed by the polytropic as well as the isothermal equation of state will be 
studied in the Newtonian, the post-Newtonian pseudo-Schwarzschild and in complete general relativistic framework. 



2 FIRST INTEGRAL OF MOTION AS A POLYNOMIAL IN CRITICAL RADIUS 



In this section, spherical and axisymmetric flows will be revisited using the Newtonian, post-Newtonian pseudo-Schwarzschild and complete 
general relativistic framework, to construct certain first integrals of motion in the form of a polynomial, a subset of real positive roots of 
which will provide the critical points of the flow. Both polytropic (for Newtonian, pseudo-Schwarzschild and general relativistic accretion in 
spherical symmetry and axisymmetry) as well as isothermal (for Newtonian and pseudo-Schwarzschild accretion in spherical symmetry and 
axisymmetry) will b e cons i dered . Some part of the c ontent of th i s section is essen tially a conglomerated overview of iDas & Sarkan d200lh : 
Sarkar & Dal j200lb ; |Pasl j2002h ; lDas et al.1 j2003l) : lDasl j2004l) : |Pas et alj J2OO7I) . We believe that such repetition will be useful to have a 



proper understanding of the subsequent sections. 
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2.1 Newtonian Accretion 

2.1.1 Spherically Symmetric Accretion 

To researches in astrophysics, physical models manifesting spherical symmetry caries an abiding appeal - studies in black hole accretion are 
no exception to this trend, since it is indeed a worthwhile exercise to consider such simplified geometric configuration in purely Newtonian 
formalism to begin with, and then to delve into more intricate problems. The pioneering work in this field was due to Bondi (1952), where 
the formal fluid dynamical equations in the Newtonian construct of space and time was introduced to study the stationary accretion problem. 

Accretion flow described in this section is 6 and <f> symmetric and possesses only radial inflow velocity. In this section, we use the 
gravitational radius r 9 as r 9 = 2GMbh /(? . The radial distances and velocities are scaled in units of r g and c respectively and all other 
derived quantities are scaled accordingly; G = c = Mbh — 1 is used hereafter. We assume the dynamical in-fall time scale to be short 
compared with any dissipation time scale during the accretion process. 

The non-relativistic equation of motion for spherically accreting matter in a gravitational potential denoted by <3? may be written as 

du du 1 dp 9$ 

— + u— + - 7 f + 7r = 0, (1) 

at or p or or 

The first term in ([TJ is the Eulerian time derivative of the dynamical velocity, the second term is the 'advective' term, the third term is the 
momentum deposition due to the pressure gradient and the last term is the gravitational force. Another equation necessary to describe the 
motion of the fluid is the continuity equation 

To integrate the above set of equations, one also needs the equation of state that specifies the intrinsic properties of the fluid. We will study 
accretion described by either a polytropic or an isothermal equation of state. 

We employ a polytropic equation of state of the form p — Kp 1 . The sound speed c s is defined by 

2 _ dp 



dp 



7-, (3) 
P 



constant entropy 

Assuming stationarity of the flow, we find the following conservation equations: 

1) Conservation of energy implies constancy of the specific energy £ 

£ = IT + + *■ (4) 

2 7 — 1 

2) Conservation of the baryon number implies constancy of the accretion rate M 

M = Am pur 1 (5) 
The corresponding entropy accretion rate comes out to be: 

2 

S = AnCs 1 ur 1 (6) 

Equation $4^ is obtained from Q}, and ([5} follows directly from (O. 

Substituting p in terms of c s and differentiating (|5j with respect to r, we obtain 

c' s = C -^pl( U - + 2 -), (7) 
2 \u r J 

where ' denotes the derivative with respect to r, Next we differentiate 10 and eliminating c' s with the help of we obtain 

u = . (8) 

u — c%/u 

A real physical transonic flow must be smooth everywhere, except possibly at a shock. Hence, if the denominator of {8} vanishes at a point, 
the numerator must also v anish at that point to e nsure the physical continuity of the flow. Borrowing the terminology from the dynamical 
systems theory (see, e.g., Jordan & Smithl 1 1999b ). one therefore arrives at the critical point conditions by making the numerator and the 



denominator of ([Sj simultaneously equal to zero. One thus finds the critical point conditions as 

(9) 









^r— r c — C sr — rc — ^ 


2 






r=r c 



The critical point and the sonic points are thus equivalent. Hereafter, the subscript c indicates that a particular quantity is evaluated at r c . The 
location of the critical point/acoustic horizon is obtained by solving the algebraic equation 

£ -U^ ± i) r ^-^= o - (io) 
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Using the system of units as described above, we substitute the value of $ in purely Newtonian construct as — 1/r and obtain: 

4£(7-l)r-(5-37)=0 (11) 

The above equation is a polynomial of degree one in r c with constant real co-efficient parametrized by [£,7]. 

We now set the appropriate limits on £ and 7 to model the realistic situations encountered in astrophysics. For relativistic flow, £ is 
scaled in terms of the rest mass energy and includes the rest mass energy, £ < 1 corresponds to the negative energy accretion state where 
radiative extraction of rest mass energy from the fluid is required. For such extraction to be made possible, the accreting fluid has to possess 
viscosity or other dissipative mechanisms, on which we would not like to focus in this work. On the other hand, although almost any £ > 1 
is mathematically allowed, large values of £ represents flows starting from infinity with extremely high thermal energy (see section 13.4 for 
further detail), and £ > 2 accretion represents enormously hot flow configurations at very large distance from the black hole, which are not 
properly conceivable in realistic astrophysical situations. Hence one sets 1 £ 2 for relativistic flows, and £ 1 for non relativistic 
flows. Now, 7 = 1 corresponds to isothermal accretion where accreting fluid remains optically thin. This is the physical lower limit for 7, 
and 7 < 1 is not realistic in accretion astrophysics. On the other hand, 7 > 2 is possible only for superdense matter with substantially 
large magnetic field (which requires the accreting material to be governed by general relativistic magneto-hydrodynamic equations, dealing 
with which is beyond the scope of this article) and direction dependent anisotropic pressure. One thus se ts 1 7 £ 2 as we ll. However, one 



should note that the most preferred values of 7 for realistic black hole accretion ranges from 4/3 to 5/3 iFrank et alj (2002). 

Eq. dl It being a polynomial of degree one, multi-criticality is ruled out in this situation. Also to note that only the real positive solutions 
of dl It in the domain r > 1 are physically acceptable, since r — 1 is the radius of the event horizon. We will observe in subsequent sections 
that due to the choice of the modified gravitational potential, or due to the departure from the spherical symmetry, the defining polynomials 
may be of degree greater than one, and the analytical solution will not be straight forward. 

For our purpose, d 1 1 b is the fundamental equation for the polytropic accretion, because its roots will provide the information about the 
multi-criticality. Irrespective of the geometric configuration, and the kind of space time in which the flow is being studied, our main aim is 
to construct the above kind of energy first integral polynomial for the critical point(s). For this purpose, the fundamental requirements are 
the expressions for the £, M and Ml. Once we have all such expressions, we can arrive at the energy polynomial, and can either solve it 
analytically (if the degree n of the polynomial is less than or equal to four) to find out the roots explicitly, or we can apply a mathematical 
scheme in general (irrespective of the degree of the energy polynomial), which will enable us to find out the number of physically acceptable 
roots (real, positive and greater than the Schwarzschild radius) corresponding to the polynomial. In this way we will be able to explore the 
transonicity of any generic flow completely analytically, which has never been done before in the literature. 

For isothermal flow, equation of state of the form 

P= P = c s p, (12) 

A 4 

is employed to study the accretion, where T is the temperature, R and fi are the universal gas constant and the mean molecular weight, 
respectively. The quantity c s is the isothermal sound speed defined by 

2 dp 



dp 



= 6T, (13) 

T 



where the derivative is taken at fixed temperature and the constant Q — Kb/ (fJ,mn) with mn — rn p being the mass of the hydrogen atom. 
In our model we assume that the accreting matter is predominantly hydrogen, hence fj, ~ 1. 
Two first integrals of motion can be obtained as: 

2 

— + QTlnp + $ = C (14) 

M = A-Kpur 2 (15) 

Where C is a constant, u and M being the dynamical flow velocity and the mass accretion rate respectively. The space gradient of the 
dynamical velocity can be expressed as: 

f = Tnr (16) 

u 

The critical point condition thus comes out to be: 



c s = er (17) 



2 

Putting the value of $ = —1/r, one obtains: 

29 2 T 2 r-l = (18) 
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which is to be solved to obtain the critical point. Once again, we obtain a polynomial of degree one, and hence the multi-criticality is ruled out. 
The generalization of the above determining equation is quite straight forward - we need to establish a polynomial relation in r c parametrized 
by the flow temperature T (for spherical accretion) as well as flow angular momentum (for axisymmetric flow). 



2. 1.2 Axisymmetric Accretion 



For all categories (Newtonian, pseudo-Schwarzschild and general relativistic) of axisymmetric accretion, We assume that the disc has a 
radius-dependent local thickne ss H(r), and its central p lane coincides wi t h the equatorial plane of t h e black hole. It is a st a ndard practice in 
accret io n disc theory (se e e.g.jMatsumoto et aljjl984h: Paczvfiski 1 1987ri ; Abramowicz et alj jl988h ; Chen & Taan] jl993h; Kafatos & Yang 



1994 h : lArtemova et al.l Jl996f) ; iNaravan. Kato & Honmj Jl997m lwiital jl998h ; lHawlev & Kroliki J200lh : lArmitage. Reynolds & Chiang 



200 lh ) to use the vertically averaged model in describing the black-hole accretion discs where the equations of motion apply to the equato- 



rial plane of the black hole. We follow the same procedure here. The thermodynamic flow variables are averaged over the disc height, i.e., a 



thermodynamic quantity y used in our model may be vertically averaged over the disc height as y = 



f" {r, (ydh) 



Jo" dh 

H(r) of the disc can be obtained by balancing the gravitational force by pressure gradient. 

For a non- viscous flow obeying the polytropic equation of state p = Kp 1 , integration of radial momentum equation: 



The local half-thickness, 



^du 1 dP d f^eff, s 
dr p dr dr \ 







(19) 



where the effective potential M is the su mmation of the gravitational potential and the 'centrifugal potential' due to the conserved 
specific angular momentum A (see iDas d2002h for further detail), leads to the following energy conservation equation (on the equatorial 
plane of the disc) in the steady state: 



1 2^ I 

—u H 

2 7 - 



— I T + * 

1 2r 2 



and the continuity equation: 



d 



— [uprHM] = 
dr 

can be integrated to obtain the baryon number conservation equation: 



(*')" 



M = x -uc s pr* 

7 



The entropy accretion rate S can be expressed as: 



— uc t 

7 



For a particular value of [£ , A, 7], it is now quite straight-forward to derive the space gradient of the acoustic velocity 
dynamical flow velocity 



dc s 
dr 



I for flow as: 



(20) 



(21) 



(22) 



(23) 



and the 



(24) 



(25) 



(26) 



rr (27) 

7 + 1 

Hence, unlike the spherical flow, the critical points and the sonic points are not equivalent. The sonic points (where the flow makes a transition 
from the subsonic to a supersonic state) are formed at a smaller (compared to the critical points) radial distance for the accretion. We will see 



dr 



7" 



and, 



du\ 
~d~r~) 



7+1 



+ $ 



dr j 
2.W 



7+1 



3_ 

2r 



1 du\ 
u dr J 



<S> ftp 



#'(r) 



2c| 
"(7+1) 



where $" represents the derivative of Hence the critical point condition comes out to be: 

$ ( r ) + 7 $ ( r ) / A 2 + r 3 $' (r) 



1+7 



3$' (r) + r$" (r 



Note that the Mach number M c at the critical point is not equal to unity, rather: 

M c = 
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in subsequent discussion that this is true in general for polytropic flow in the axisymmetric configuration if the disc height is taken not to be 
a constant quantity. 

For any fixed set of [£ , X, 7], the critical points can be obtained by solving the following polynomial of r: 



£ - 



2r 2 



27 



$ ( r ) + 7 $ ( r ) / A 2 + r 3 $ (r) 



3$'(r) + r$"(r) 



= 0. 



(28) 



Similarly, for the isothermal accretion, the critical point conditions comes out to be (see 
detail): 
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2003h or the derivation and further 



QT2 



\ 



2r, 



(29) 



Note that the Mach number at the critical point is exactly equal to unity, hence, unlike the polytropic axisymmetric flow, the critical points 
and the sonic points are identical for isothermal accretion disc. 

The solution of the following equation parametrized by [T, A] provides the critical/sonic points: 



2 2 

Q2y )r=r 



+ 



2A 2 



r s e TQ 2 r 3 h 



= 



(30) 



Putting the Newtonian potential (4> = — 1/r) in (28) and l |30t , we find that for polytropic as well as for the isothermal flow, the deter- 
mining polynomials are quadratic equations: 



2£r c + 2 



2O7 
~ 1 



r c + A 2 



2O7 
~ 1 



Polytropic 



50 2 Tr 2 + 2r c - 2A 2 = Isothermal 



(3D 
(32) 



Hence, the maximum number of the critical points such a flow can have is two. However, real physical flow does not allow such configuration. 
The theory of dynamical systems asserts that no two successive critical points can be of similar type (both saddle, or both centre type, for 
exam ple). Hence, out of the two solut i ons, flow will have one saddle type and one centre type critical point. H owever, it has been demonstrated 
that I Muc hotrzeb & Paczvf iski 1982; Muchotrzeb 1983: M uchotrzeb-Czern v 1986: Chaudh urv et alj200a) a saddle type critical point which 
is accompanied by a centre-type critical point but without another saddle type critical point, the flow solutions will all curl about the centre- 
type critical point by forming a homoclinic orbit passing through the saddle type critical point, and by no means connecting the event horizon 
with infinity through such solutions will be possible. To avoid this globally invalid situation, and to make accretion a feasible proposition, 
purely Newtonian axisymmetric flow must be mono-transonic. This finding also enable us to introduces a generic statement of the form: 

Proposition 2.1. No physically acceptable accretion flow can have two critical points. 

Hereafter in this section we provide (after citing proper references) only the expressions for the critical point conditions and the form of the 
corresponding first integral polynomials for the critical point(s), parametrized by [£,y] (spherical flow) or [£, A, 7] (axisymmetric flow) for 
polytropic accretion, and by [T] (spherical flow) or [T, A] (axisymmetric flow) for the isothermal accretion. 



2.2 Post-Newtonian pseudo-Schwarzschild Accretion 



Rigorous investigation of the complete general relativistic transonic black hole accretion disc structure is extremely complicated. At the same 
time it is understood that, as relativistic effects play an important role in the regions close to the accreting black hole (where most of the 
gravitational potential energy is released), purely Newtonian gravitational potential cannot be a realistic choice to describe transonic black 
hole accretion in general. To compromise between the ease of handling of a Newtonian description of gravity and the realistic situations 
described by complicated general relativistic calculations, a series of 'modified' Newtonian potentials have been introduced to describe the 
genera l relativistic effects that are most important for accretion disk structure around Schwarzschild and Kerr black holes (see lArtemova et alj 
1996h:lDasli [2002J), and references therein for further discussion). 



Introduction of such potentials allows one to investigate the complicated physical processes taking place in disc accretion in a semi- 
Newtonian framework by avoiding pure general relativistic calculations so that most of the features of spacetime around a compact object are 
retained and some crucial properties of the analogous relativistic solutions of disc structure could be reproduced with high accuracy. Hence, 
those poten tials describing the spa c e time arou nd a non rotating black hole might be designated as 'post-Newtonian pseudo- Schwarzschild' 



potentials' (lArtemova et al] dl99a) : ] 

close (within 1 — 2 gravitational radius) to a black hole horizon, it is safe to use the pseud o-Schwarzschild potential s to study accretion on 



Das (2002), and references therein). As long as one is not interested in astrophysical processes extremely 



to a Schwarzschild black hole. Among all available pseudo-Schwarzschild potentials, the 1 Paczvnski & Wiitj|l980l) ) potential of the form 
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-1/2 (r — 1) serves the best for the aforesaid purpose JArtemova et al.|[T996l Das & Sarkaj2001 : Das[ 2002). and in this paper we will 



study the transonic structure under the influence of this potential only. 
2.2.1 Spherical Accretion 

Following the procedure described in 32.1.11 the corresponding polynomials in r c parametrized by [£,j] and [T] for the polytropic and 
isothermal flow respectively can be expressed as: 

8£( 7 - l)rl + (16£ - 16£"7 + 37- 5)r c + (8^7- S£ - 4 7 + 4) = (Polytropic) (33) 
4r 2 9 2 T 2 - (1 + 89 2 T 2 )r c + 46 2 T 2 = (Isothermal) (34) 



According to the proposition | |2.U . the only physically acceptable flow will be mono-transonic for spherical accretion under the influence of 
the jPaczvnski & WiitJl980h potential. 

2.2.2 Axisymmetric Accretion 

Following the procedure described in 32.1.21 we find that for the polytropic accretion the energy first integral polynomial is of the form: 

a4,rt + a 3 rl + a 2 r 2 + a\r c + do = (35) 
where: 

a.4 = -10£ ,a 3 = 16£ - 5 + 2 ^ - , a 2 = -6£ + 3 + 5A 2 - 4 T> ,ai = y———rr , 0,0 = -A 2 ( 7 + | ) (36) 

(7-1) (7-I) (7-1) 7-1 

The first integral polynomial for isothermal flow is: 

2r* - (2 + 9 2 T)r^ + (39 2 T)r 2 - 2A 2 r c + 2A 2 = (37) 

Both of the above polynomials are of degree n = 4. However, we will prove in the next section and in Sj4]that at most three physically 
acceptable (real, positive and greater than unity) roots may be found. Hence the accretion flow governed by i35l or d37b is multi-critical. 

2.3 Analytical solutions 

Subjected to a set of specific initial boundary conditions - a set of values of [£, A, 7] for the adiabatic flow and a set of values of [T, A] 
for the isothermal flow, equation d35t and J37b can be solved completely analytically to find the number of real physical roots. We use the 
Ferrari's method for solving a quartic equation to find out the solutions for ([35} and J37b . The Ferrari's method has been described in detail 
in Appendix - 1. 

As for an illustrative example, a choice of [£ — 0.005, A = 1.65, 7 = 4/3] provides four roots for the equation d36t as 51.54, 6.3713, 
2.8315 and 0.75, all measured in the unit of 2GM/c 2 . The fourth root is located inside the event horizon, hence the accretion flow can have 
three real physical critical points. If one uses the eigenvalue analysis I Chaudhurv et al J 20061 : Goswami et alj2007 ) to find out the nature of 



the critical points, then it can easily be shown that the outer and the inner saddle type critical points are located at 51.54 and 2.8315 and the 
centre type critical point is located at 6.3713. For isothermal flow, [T = 2x 10 9 , A = 1.9] provides four roots as 205.198875, 4.328, 2.2689 
and 0.7566 (all measured in the unit of 2GAl/c 2 ). Here again, the fourth root is located inside the event horizon, hence the accretion flow 
can have three real physical critical points. The above mentioned eigenvalue analysis can reveal the nature of the critical points, and once 
again it can be shown that there are two saddle type inner and the outer critical points and one centre type middle critical point. As mentioned 
earlier, for the isothermal flow the sonic points are identical with the critical points, but for the polytropic flow they are not, and one needs to 
numerically integrate the flow from the saddle type critical points to find the corresponding sonic points (the radial distance where the Mach 
number becomes exactly equal to unity). 

2.4 General Relativistic flow in the Schwarzschild Metric 

2.4.1 Spherical Flow 

The critical point condition for such flow comes out to be 1 Dasll2004 ): 



4r c -3 

Thus the critical points and the sonic points are isomorphic. The corresponding energy first integral polynomial is: 



(38) 



A 3 r 3 + A 2 r 2 +Air + A Q = (39) 



Transonicity using Sturm chains 9 



where 

A = 64(7 - l) 2 - 64£V - 64£ 2 + 1287^ 2 



At 

As 



144(7 
188(7 - 
-27(7 - I) 2 



l) 2 + 196£ 2 7 2 + 64£ 2 - 160£ 2 7 



1) 



36fV 



16£ 2 



48£ 2 7 



(40) 



However, it can be shown that inspite of the fact that the above equation is a cubic one, general relativistic spherical accretion is only 
mono-transonic (Das 2004). 



2.4.2 Axisymmetric Accretion 

Using r g = GMbh/c 2 , the critical point conditions are jPas et alj2007l) 



u c = ± 



/ 2 (r c ,A) _ 7 +l 



/i(rc,A) + /2(r c ,A)' 



h(ro,X) 



A(rc,A) 



where u c = u(r c ) and c c = c s (r c ), r c being the location of the critical point. fi(r c , A) and /a(r c , A) are defined as: 



(41) 



, 3rg - 2A 2 r e + 3A 2 2r c - 3 2r e 3 - A 2 r c + A 2 

/l(rc ' A) = ^-A»r e (r e -2)' = 2) " J? - AV c (r c - 2) (42) 

Clearly, the critical points are not coincident with the sonic points. 
We derive the energy first integral polynomial as: 

rl 4 {(l - 1) 2 (-108)} + r c 13 {396( 7 - l) 2 } 
+r c 12 {( 7 -l) 2 (252A 2 - 360)} 
+^ 1 {(7 - 1) 2 (-1356A 2 )} 
+r c 10 {(2424)( 7 - 1) 2 A 2 + 3£ 2 ( 7 + l) 2 
+432£ 2 ( 7 - l) 2 - 192(7 - 1) 2 A 4 } 
+^c{(7 - 1) 2 (1352A 4 - 1440A 2 ) - £ 2 (j + 1) 2 (6A 2 + 6)36£ 2 (7 2 - 1) - 864£ 2 (7 - l) 2 } 
+^{(7 ~ 1) 2 (48A 6 - 3598A 4 ) + £ 2 (j + 1) 2 (3A 4 + 34A 2 ) + 36£ 2 (7 2 - 1)A 2 - 864A 2 £ 2 (7 - l) 2 )} 
+rl{( 7 - 1) 2 (4160A 4 - 416A 6 ) - £ 2 ( 7 + 1) 2 (65A 2 + 26A 4 ) - 96£ 2 (7 2 - 1)A 2 + 3024£ 2 ( 7 - 1) 2 A 2 } 
+r 6 c {{ 1 - 1) 2 (1448A 6 - 1800A 4 ) + £ 2 ( 7 + 1) 2 (90A 4 + 42A 2 ) - £ 2 ( 7 2 - 1)(48A 4 - 72A 2 ) + £ 2 ( 7 - 1) 2 (576A 4 - 2592A 2 )} 
+r 5 c {(j - 1) 2 (-2512A 6 ) + £ 2 (j + 1) 2 (23A 6 - 156A 4 ) + 248£ 2 ( 7 2 - 1)A 4 - 2880£ 2 A 4 ( 7 - l) 2 } 
+r 4 {2152A 6 ( 7 - l) 2 + £ 2 ( 7 + 1) 2 (136A 4 - 30A 6 ) - £ 2 ( 7 2 - 1)(432A 4 - 16A 6 ) + £ 2 ( 7 - 1) 2 (4800A 4 - 128A 6 )} 
+r 3 {( 7 - 1) 2 (-720A 6 ) + £ 2 ( 7 + 1) 2 (280A 6 - 48A 4 ) - £ 2 ( 7 2 - 1)(112A 6 - 252A 4 ) + £ 2 ( 7 - 1) 2 (832A B - 2592A 4 )} 

+r 2 {£ 2 ( 7 + 1) 2 (-400A 6 ) + 292£ 2 A 6 ( 7 2 - 1) - 2048£ 2 ( 7 - 1) 2 A 6 } 
+r c {304e 2 A 6 ( 7 + l) 2 - 336£ 2 A 6 (7 2 - 1) 
+2208£ 2 A 6 ( 7 - l) 2 } + {-96£ 2 A 6 ( 7 + l) 2 + 144£ 2 A 6 ( 7 2 - 1) - 864£ 2 A 6 ( 7 - l) 2 } = (43) 

The above equation, being an n = 14 polynomial, is non analytically solvable. In Sj4]we will demonstrate how we can analytically find out 
the number of physically admissible real roots for this polynomial, and can investigate the transonicity of the flow. 



3 STURM THEOREM AND GENERALIZED STURM SEQUENCE (CHAIN) 

In this section we will elaborate the idea of the generalized Strum sequence/chain, and will discuss its application in finding the number of 
roots of a algebraic polynomial equations with real co-efficients. Since the central concept of this theorem is heavily based on the idea of the 
greatest common divisor of a polynomial and related Euclidean algorithm, we start our discussion by clarifying such concept in somewhat 
great detail for the convenience of the reader. 



3.1 Greatest common divisor for two numbers 

Let us start with the concept of the divisibility first. 

Given two non-zero integers Z\ and za, one defines that z\ divides 22, if and only if there exists some integer 23G Z such that: 
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Z2 — Z3Z2 (44) 
The standard notation for the divisibility is as follows: 

Zi|z2means 'zi divides z' 2 (45) 

The concept of divisibility applies to the polynomials as well, we treat such situations in the subsequent paragraphs. 

Now consider two given integers z\ and Z2, with at least one of them being a non-zero number. The 'greatest common divisor' (or the 
'greatest common factor' or the 'highest common factor') of Z\ and Z2, denoted by gcd(zi, Z2), is the positive integer Zd£^, which satisfies: 

i)zd\zi and Zd\z\. 
ii)For any other z c £Z, if z c \z\ andz c |z2 

thenz c jz d (46) 

In other words, the greatest common divisor gcd(zi, 22) of two non zero integers zi and 22 is the largest possible integer that divides both 
the integers without leaving any remainder. Two numbers z\ and Z2 are called 'co-prime] (alternatively, 'relatively prime'), if: 

gcd(z 1 ,z 2 ) = l (47) 

The idea of a greatest common divisor can be generalized by defining the greater common divisor of a non empty set of integers. If Sz is a 
non-empty set of integers, then the greatest common divisor of Sz is a positive integer Zd such that: 

i) If Zd|zifor all zi£Sz 

ii)If Z2I21, for all zi£Sz, then Z2\zd (48) 
then we denote Zd ~ gcd(Sz). 



3.2 Euclidean algorithm 

Euclidean algorithm (first described in detail in Euclid's 'Elements' in 300 BC, and is still in use, making it the oldest available numerical 
algorithm s till in common use) provides an efficient procedure for computing the greatest common divisor of two integers. Following Stark 
i Starklll978h . below we provide a simplified illustration of the Euclidean algorithm for two integers: 



Let us first set a 'counter' i for counting the steps of the algorithm, with initial step corresponding to i = 0. Let any ith step of 
the algorithm begins with two non-negative remainders r%—i and ri-2 with the requirement that n—i < n-2, owing to the fact that the 
fundamental aim of the algorithm is to reduce the remainder in successive steps, to finally bring it down to the zero in the ultimate step which 
terminates the algorithm. Hence, for the dummy index i, at the first step we have: 

r_ 2 = z 2 and r_i = z\ (49) 

the integers for which the greatest common divisor is sought for. After we divide Z2 by z\ (operation corresponds to i = 1), since Z2 is not 
divisible by zi, one obtains: 

r-2 = qor-i + ro (50) 

where ro is the remainder and go be the quotient. 

For any arbitrary ith step of the algorithm, the aim is to find a quotient qj and remainder r», such that: 

n-2 — qtn-i +ri, where r; < ri-i (51) 

at some step i = j (common sense dictates that j can not be infinitely large), the algorithm terminates because the remainder becomes zero. 
Hence the final non-zero remainder tj-i will be the greatest common divisor of the corresponding integers. 

We will now illustrate the Euclidean algorithm for finding the greatest common divisor for two polynomials. 



3.3 Greatest common divisor and related Euclidean algorithm for polynomials 

Let us first define a polynomial to be 'monic' if the co-efficient of the term for the highest degree variable in the polynomial is unity (one). 
Let us now consider pi(x) and P2{x) to be two nonzero polynomials with co-efficient from a field F (field of real, complex, or rational 
numbers, for example). A greatest common divisor of pi (x) and p2 (x) is defined to the the monic polynomial Pd(x) of highest degree such 
that pd(x) divides both pi(x) and P2(x). It is obvious that F be field and Pd(x) be a monic, are necessary hypothesis. 

In more compact form, a greatest common divisor of two polynomials pi , p 2 £M,[X\ is a polynomial pd£R[X] of greatest possible degree 
which divides both pi and p 2 . Clearly, pd is not unique, and is only defined upto multiplication by a non zero scalar, since for a non zero scalar 
cGR, if pd is a gcd(pi,p2&R\X]), so as cpd- Given polynomials pi,p2GK[X], the division algorithm provides polynomials p3,p4£R[X], 
with deg(p<i) < deg(ps) such that 
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Pi = P3P2 + Pi (52) 

Then, if pd is gcd(pi,p2), if and only if pd is gcd(p2,PA) as is obvious. 

One can compute the gcd of two polynomials by collecting the common factors by factorizing the polynomials. However, this technique, 
although intuitively simple, almost always create a serious practical threat while making attempt to factorize the large high degree polynomials 
in reality. Euclidean algorithm appears to be relatively less complicated and a faster method for all practical purposes. Just like the integers 
as shown in the previous subsection, Euclid;s algorithm can directly be applied for the polynomials as well, with decreasing degree for the 
polynomials at each step. The last non-zero remainder, after made monic if necessary, comes out to be the greatest common divisor of the 
two polynomials under consideration. 

Being equipped with the concept of the divisibility, gcd and the Euclidean algorithm, we are now in a position to define the Strum 
theorem and to discuss its applications. 



3.4 Sturm Theorem: The purpose and the definition 

Sturm theorem is due to Jacaues Charles Francois Strum {29 th September 1803 - 15th December 1855), a Geneva born French mathematician 
and a close collaborator of Joseph Liouville (Strum was the co-eponym of the great Sturm-Liouville problem, an eigenvalue problem in sec- 
ond order differential equation). In collaboration with his long term friend Jean-Daniel Colladon, Sturm performed the first ever experimental 
determination of the velocity of sound in water. Sturm theorem, published in 1829 in the eleventh volume of the 'Buletin des Sciences de 
Ferussac' under the title 'Memoire sur la resolution des equations numeriquesQ. Sturm theorem, which is actually a root counting theorem, 
is used to find the number of real roots over a certain interval of a algebraic polynomial with real co-efficient. It can be stated as: 

Theorem 3.1. The number of real roots of an algebraic polynomial with real coefficient whose roots are simple over an interval, the endpoints 
of which are not roots, is equal to the difference between the number of sign changes of the Sturm chains formed for the interval ends. 

Hence, given a polynomial p€R[X], if we need to find the number of roots it can have in a certain open interval ]a,b[, a and b not being 
the roots of /, we then construct a sequence, called 'Sturm chain', of polynomials, called the generalized strum chains. Such a sequence is 
derived from p using the Euclidean algorithm. For the polynomial p as described above, the Sturm chain po , p\ . . . can be defined as: 

Po = P 
Pi = P 

p n = -rem(p„-2,Pn-i) ,n>2 (53) 

where rem (p n -2,Pn-i) is the remainder of the polynomial p„-2 upon division by the polynomial p n -i- The sequence terminates once 
one of the pi becomes zero. We then evaluate this chain of polynomials at the end points a and b of the open interval. The number of roots 
of p in ]a, b[ is the difference between the number of sign changes on the chain of polynomials at the end point a and the number of sign 
changes at the end point b. Thus, for any number t, if N p (t) denotes the number of sign changes in the Sturm chain po(t), pi (t), then for 
real numbers a and b that (both) are not roots of p, the number of distinct real roots of p in the open interval ]a, b[ is [-W p ( a ) — -Wp(ft)] • By 
making a^— oo and 6— >+oo, one can find the total number of roots p can have on the entire domain of K. 
A more formal definition of the Strum theorem, as a corollary of the Sylvester's theorem, is what follows: 

Definition Let R be the real closed field, and let p and P be in R[X].7he Sturm sequence of p and P is the sequence of polynomials 

(po,pi, ...,Pfe) defined as follows: 

Po =p,pi= p'P 

Pi = Pi-iQi — fi-2 with qi £ R[X] and deg(p.;) < deg(pi-i) for i = 2, 3, k, Pk is a greatest common divisor of p and p'P. 

Given a sequence (an, a-k) of elements of R with an 7^ 0, we define the number of sign changes in the sequence (ao, ■ ■■ak) as 
follows: count one sign change if a^a; < with I > i + 1 and dj = for every j, i < j < I. 

If a G R is not a root of p and (po, Pk) is the Sturm sequence of p and P, we define v(p, P; a) to be the number of sign changes in 
(po(a),...p k (a)). 

Theorem 3.2. (Sylvester's Theorem) Let R be a real closed field and let p and P be two polynomials in R[X]. Let a, b £ Rbe such that 
a < b and neither a nor b are roots of p. Then the difference between the number of roots of p in the interval ]a, b[ for which P is positive 
and the number of roots of p in the interval ]a, b[ for which P is negative, is equal to v(p, P; a) — v(p, P; b) 

Corollary 3.3. (Sturm's Theorem): Let R be a real closed field and p £ R[X\. Let a,b £ R be such that a < b and neither a nor b are 
roots of p. Then the number of roots of p in the interval ]a, b[ is equal to v(p, 1; a) — v(p, 1; b). 

The proof of these two theorems are given in the Appendix. 

In next section we will provide an illustrative application of Sturm theorem by applying it to compute the number of critical points for 
axisymmetric polytropic accretion in Paczynski and Wiita (1980) pseudo-Schwarzschild potential. 



1 According to some historian, the theorem was originally discovered by jean Baptist Fourier, well before Sturm, on the eve of the French revolution. 
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4 NUMBER OF CRITICAL POINTS FOR PSEUDO-SCHWARZSCHILD POLYTROPIC AXISYMMETRIC ACCRETION 

In i]2.2.1l we have shown that the energy first integral of motion polynomial in critical points is a fourth degree equation for this case, and 
hence analytical solution is possible to explicitly find the roots (which has been performed in ^2.3b . We compute the number of roots (for 
same set of initial boundary condition) as a sanity check against the known (already calculated, both analytically as well as numerically) 
to show that Sturm theorem can actually provide the exact number of sonic point, before we apply it to find out the number of roots for 
relativistic accretion where no such analytical solution exists (since it is a fourteenth degree equation as already been demonstrated in i)2.4.2b 
for cross verification. 

The energy first integral polynomial for such situation can be expressed as: 

po(r) = a 4 r 4 + a 3 r 3 + a 2 r 2 + a±r + a (54) 

where 

£14 = — 10£ , 

a 3 = 16£-5 + 27 (55) 
(7-1) 



a 2 = -6£ + 3 + 5A - 



ai 



2 4A 2 7 



(7-1)' 

8A 2 
(7-1)' 



ao = -A 2 (^— ) (56) 
7-1 

Subsequently we obtain the corresponding sequences as: 

pi(r) = 4d4r 3 + 3a 3 r 2 + 2a 2 r + ai 
Pa(r) = hr 2 + b 2 r + b z 
Pz(r) = cir + c 2 

P4(r) = -[6 3 -(6 2 -— )-] (57) 

Cl Ci 

where 



r02 


3(a 3 ) 2 




16a4 


3ai 


2a 2 a 3 




16a4 



, r aia 3 

<4= _Pa a _^_ (3aB _^)^] > 
oi oi 6i 

r , 4a 4 6 2 . & 3 , 
c 2 = - a a - (3a 3 — ) — 

01 Ol 



For the parameter set 



£ = 0.0001, A = 1.75,7= - 



, one obtains 



Po(r 
Pi{r 
p 2 (r 
Ps{r 
Pi(r 



-O.OOlr 4 + 3.00016r 3 - 30.6881r 2 + 73.5r - 39.8125 

-0.004r 3 + 9.00048r 2 - 61.3762r + 73.5 

-1672. 256r 2 + 11453.475r - 13742.088 

-0.117r + 0.2388 

-2668.4868 



We now consider the left boundary of the open interval to be the Schwarzschild radius, which is unity in our scaled unit of radial distance 
used to describe the flow equation for pseudo-Schwarzschild axisymmetric polytropic accretion, and the right boundary to be 10 6 (in our 
scaled unit of radial distance), which is such a large distance that beyond which practically no critical point is expected to form. We then 
calculate the sign changes of the Sturm sequence and show the results in the following table: 
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r = 1 


r = 


10 6 


po(l) = +ve 


po(10 6 ) = 


—ve 


pi(l) = +«e 


Pi(10 6 ) = 


—ve 


p 2 (l) = -«e 


p 2 (10 6 ) = 


—ve 


P3(l) = 


p 3 (10 6 ) = 


—ve 


p 4 (l) = -«e 


P4 (io 6 ) = 


—ve 



We find that iV/m = 3 and Nfn 6) = 0, hence the number of critical point obtained is three, which is fully consistent with previous 
analytical and numerical computation. 

In next section, we illustrate how one can apply Sturm theorem to find out the number of critical point for general relativistic axisym- 
metric accretion in Schwarzschild metric. 



5 NUMBER OF CRITICAL POINT FOR RELATIVISTIC ACCRETION 

As has been demonstrated for the pseudo-Schwarzschild case, we first write down the complete expression for the Sturm chains. Then for a 
suitable parameter set [£, A, 7], we can find the difference of the sign change of the Sturm chains at the open interval left boundary, i.e., at 
the Schwarzschild radius, and at the right boundary, i.e., at 10 6 gravitational radius, to find the number of critical points the accretion flow 
can have. 

The form of the original polynomial has already been explicitly expressed using d43t . We now construct the Sturm chains as: 

po(r) = A 14 r 14 + A w r 13 + A 12 r 12 + ....+ A x r + A 

pi(r) = UA 14 r ri + 13A ls r 12 + .... + 2A 2 r + Ax 

P2(r) = Bi3r 12 + B 12 r 11 + ....+ B 2 r + Bi 

Ps(r) = ^r 11 +C lir 10 + .... + C 2 r + Cx 

p 4 (r) = D lir 10 +D 10 r 9 + .... + D 2 r + Dx 

Ps(r) = E 10 r 9 + E 9 r s + .... + E 2 r + Ex 

p 6 (r) = F 9 r s + F a r 7 + ....+ F 2 r + F 1 

p 7 (r) = G»r 7 + G 7 r 6 + .... + G 2 r + Gx 

ps(r) = H 7 r 6 + H 6 r 5 + ....+H 2 r + H 1 

P9{r) = Ier 5 + I 5 r 4 + ....+ I 2 r + h 

Pxo(r) = J 5 r 4 + Jxr 3 + .... + J 2 r + Jx 

pxx (r) = K 4 r 3 + K 3 r 2 + K 2 r + Kx 

Px 2 {r) = L 3 r 2 + L 2 r + Lx 

pi 3 (r) = M 2 r + Mx 

n rr IT L3M1 Mi, 

Pl4 (r) = - { Lx-(L 2 -—) — } 

Where the explicit expression of the corresponding co-efficients Ai,Bi... has been provided in the equation ([43} and in the Appendix - III. 

It is important to note that direct application of the Sturm's theorem may not always be sufficient since some of the roots may yield a 
negative energy for £ (since the £ equation was squared to get the polynomial). To get positive values of the energy, we must impose the 
condition that 

y-(l + c 2 )>0, (58) 
which is the term present in £ which could go negative. This introduces the condition that 

2 7r 4 - 2r 4 - 5r 3 7 + 3r 3 - 3AV 7 + 5A 2 r 2 + 10A 2 r7 - 18A 2 r - 8A 2 7 + 16A 2 



2(r 4 - 2A 2 r 2 + 7A 2 r - 2r 3 - 6A 2 ) 



< (59) 



p(v) 

This is equivalent to the condition that ; ^ where p(r) and q[f) are 4th order polynomials. To find the region where this happens, one 

q{r) 

has to find the 4 roots of each of p(r) and q(r) - which is analytically possible since roots of quartics are analytically solvable. Once the 
roots are obtained it is a trivial matter to check for what regions the rational function is positive. 

A simplified version for the above mentioned procedure to find the positivity condition is as follows: 
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We would like to find out the intervals in which p(r)/q(r) > where p(r) and g(r) are quartic polynomials. We factorize p(r) — 
(r — ri)(r — r2)(r — rs)(r — rs) and q(r) = (r — si)(r — S2)(r — s;i)(r — S4) using the algorithm for finding roots of a quartic. If the 
roots are all real, we can just note down the sign changes of each factor to the right and left of each root and find out the intervals where the 
rational function is positive. If there are complex roots, they come in complex conjugates, since the coefficients of the polynomials are real. 
Say, if rs is complex and n is its complex conjugate, then the part (r — rs)(r — n) — r 2 — (rs + n)r + rsn does not change sign since 
it is non-zero on the real line. It is easy to determine its sign. 

As an illustrative example we take [£ — 1.0001, A = 3.4, 7 = 4/3]. We then get the number of sign changes to be three in the interval 
between the location of the event horizon and a representative distance 10 6 (in the units of GM/c 2 ). Hence the accretion is multi- critical for 



this s et of initial boundary condition. This is in complete agreement with the results found elsewhere (see, e.g., the figure 2. of dDas et al 



2007)) using numerical techniques. 



6 DISCUSSION 

Our methodology is based on the algebraic form of the first integral obtained by solving the radial momentum equation (the Euler equation 
to be more specific, since we are confined to the inviscid flow only). The structure for su ch a first integral has to be a formal polynomial 
with appropriate constant co-efficients. There are certain pseudo-potentials, introduced by I Artemova et al. IHiil), for example, for which 



the first integrals (neither for the polytropic nor for the isothermal flow) can not be written in a purely algebraic polynomial form. Also for 
general relativistic accretion in the Kerr metric, the expression for the energy first integral can not be reduced to such a polynomial form. 
Hence, the Sturm's generalized chain can not be constructed for such accretion flow. Alternative methodology are required to investigate the 
multi-critical behaviour for such kind of accretion. 

Using the method illustrated in this work, it is possible to find out how many critical points a transonic black hole accretion flow can 
have. It is thus possible to predict whether such accretion flow can have multi-critical properties for a certain speci fic value/domain of the 
initial boundary cond itions. It is, however, not possible to investigate, using the eigenvalue analysis as illustrated in Jchaudhurv et alJbood : 



dary 



Goswami et alj |2007). the nature of such critical points - i.e., whether they are of saddle type or are of centre type, since such prediction 
requires the exact location of the critical points (the value of the roots of the polynomial). However, the theory of dynamical systems ensures 
that no two consecutive critical points be of same nature (both saddle or both centre). On the other hand, our experience predicts (it is 
rather a documented fact) that for all kind of black hole accretion, irrespective of the equation of state, the space time geometry or the flow 
configuration used, one has two saddle type critical points and one centre type critical point flanked by them. Hence if the application of 
Sturm's generalized chain ensures the presence of three critical points, we can say that out of those three critical points, accretion flow will 
have two saddle type critical points, hence a specific subset of the solution having three roots corresponding to the first integral polynomial, 
can make transonic transition for more than one times, if appropriate conditio ns for c onnecting the flow through the outer critical point and 



can make transonic transition tor more man one times, it appro priate condition s tor connecting the now tnr 
for flow through the inner critical points are available, see, e.g jDas & Czernv for further discussion 



In this work we have considered only inviscid accretion. Our methodology of investigating the multi-critical properties, however, is 
expected to be equally valid for the viscous accretion disc as well. For the viscous flow, the radial momentum conservation equation involving 
the first order space derivative of the dynamical flow velocity will certainly provide a first integral of motion upon integration. Because of the 
fact that a viscous accretion disc is not a non-dissipative system, such constant of motion, however, can never be identified with the specific 
energy of the flow. The integral solution of the radial momentum equation would then be an algebraic expression of various flow variables 
and would perhaps involve certain initial boundary conditions as well. Such an algebraic expression would actually be a constant of motion. 
What exactly would that expression physically signify, would definitely be hard to realize. However, one may perhaps arbitrarily parametrize 
that conserved algebraic expression using some astrophysically relevant outer boundary conditions, and if such algebraic expressions can 
finally be reduced, using the appropriate critical point conditions, to an algebraic polynomial form of the critical points, construction of a 
generalized Sturm chain can be made possible to find out how many critical points such an accretion flow can have subjected to the specific 
initial boundary condition. Since for accretion onto astrophysical black holes, having multiple critical points is a necessary (but not sufficient) 
condition to undergo shock transition, one can thus analytically predict, at least to some extent, which particular class of viscous accretion 
disc are susceptible for shock formation phenomena. 

Our work, as we believe, can have a broader perspective as well, in the field of the study of dynamical systems in general. For a first order 
autonomous dynamical system, provided one can evaluate the critical point conditions, the corresponding generalized nth degree algebraic 
equation involving the position co-ordinate and one (or more) first integral of motion can be constructed. If such algebraic equation can 
finally be reduced to a nth degree polynomial with well defined domain for the constant co efficient, one can easily find out the maximal 
number of fixed points of such dynamical systems. 
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7 APPENDIX - 1 : ON FERRARI'S METHOD 

Ferrari's method is used to analytically solve a quartic equation for its roots. Given the quartic equation Ax* + Bx 3 + Cx 2 + Dx + E = 0, 
(A, B, C, D real or complex) its solution (i.e. the roots of the quartic) can be found by means of the following calculations: 

a _ B 3 _ BC , D 

> J SA 3 ' 2A? "r" A 

_ 3B 4 , CB 2 _ BP , E 

' 256A 4 " 1 ~ + A 



If/3 = then x = ± s ^ z2*a£EE , 

where ± s and± t are two distinct sets of plus and minuses , i.e., there are four possibilities, (± s , ± t ) = (+,+), (+, —),(—,+),(—,—). 
If P / then continue with 

P = -— - -Y 

^ 12 r> 

n — « 3 "7 /3 2 



R = — 2 j- ^1 Ql 4- ^ (either sign of the square root will do) 

[/ = R 1 ^ 3 (there are three complex roots, but any one of them will do). 

To compute the cube-root of a complex number, we proceed as follows: Let R = xi + 1x2 = re %6 where r = \J x\ + x\ and 
d = tan- 1 (^). (Here i = 

R 1/3 = r 1/3 e e/3 = r 1/3 (cos(8/3)+ism{9/3)). 

The other possible choices are R 1 ^ 3 ^ and R 1 ^ 3 ^ 2 where ui = (— 1 + \/3i)/2 is a cube root of unity. 
If U = let y = - Q 1/3 
UU£Q]et y = - Ss. + U-JL 
Let = V a + 2 y 
Then, the four roots of the quartic are 

X ~ 4A + 2 

The two ± s must have the same sign and the ±t are two independent ±. To get all the roots, compute x for the four possibilities of ± s 
and ±t- 



8 APPENDIX - II : PROOF OF SYLVESTER'S THEOREM: 

First note that the Sturm sequence (/o, ---fk) is (up to signs) equal to the sequence obtained from the euclidean algorithm. Define a new se- 
quence (g , ...,gk) by gi = fi/fk fori G {0, k}. Note that the number of sign changes in (fo(x), fi(x)) (resp. (/i_i(a:), fi(x), f i+1 (x))) 
and the number of sign changes in (go(x), gi{x)) (resp. (gi-i(x), gi(x), gi+i(x))) coincide for any x which is not a root of /. Note also 
that the roots of go are exactly the roots of / which are not roots of g. Observe that for i £ 0, k,gi-i and gi are relatively prime. We 
consider, now, how v(f,g; x) behaves when x passes through a root c of a polynomial gi. If c is a root of go, then it is not a root of gi. We 
write /'(c) > (resp. < 0) if /' is positive ( resp. negative ) immediately to the left of c. The sign of /'(c+) is defined similarly. Now we 
recall the following result: if R is a real closed field, / G R[X],a, b £ R with a < b and if the derivative /' is positive (resp. negative) on 
]o, 6[, then / is strictly increasing (resp. strictly decreasing) on [a, b]. Then, according to the signs of g(c), /'(c_) and /'(c+) we have the 
following 8 cases: 

S (c)>0,/'(c_)>0,/'(c+)>0 





c c + 


/ - 


+ 


f'g + 


+ 



S (c)<0,/'(c_)>0,/'(c+)>0 
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C- c c+ 

f - + 
f'g - r 



g(c)>0,f'(c-) <0,/'(c+) >0 





c c + 


/ + 


+ 


/'s - 


+ 



S(c)<0,/ , (c_) <0,/'(c+) >0 



c_ c c + 


/ 


+ + 




+ 



g(c)>0,f'(c-)>0,f'(c+)<0 



C- c c + 

f - o - 



S (c)<0,/'(c_)>0,/'(c+)<0 



c_ c c + 
/ - - 



S (c)>0,/'(c_)<0,/'(c + )<0 



c_ c c + 
/ + - 



S (c)<0,/'(c_)<0,/'(c+)<0 
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C- 


c c+ 


/ 


+ 


- 


f'g 


+ 


+ 



In every as x passes through c, the number of sign changes in (fo(x),fi (x)) decreases by 1 if g(c) > 0, and increases by 1 if g(c) < 0. 
If c is a root of g% with i = 1, ...k, then it is neither a root of nor a root of gi+i, and 3;-i(c)g;+i(c) < 0, by the definition of the 
sequence. Passing through c does not lead to any modification of the number of sign changes in (fi-i(x), fi(x), fi+i(x)) in this case. 

Proof of Sturm's theorem: Using g — 1 in previous theorem. 



9 APPENDIX - III: EXPLICIT EXPRESSIONS FOR THE CO-EFFICIENTS FOR THE STURM CHAIN CONSTRUCTED 
FOR THE RELATIVISTIC AXISYMMETRIC ACCRETION 

2A 12 13A? 3 1 

D 3An 12^ 12 A 13 
Bl2 = _{ ^4 196^ } 

Bn - " { ^4 196^r } 

Bl ° = ^IT - 196A 14 1 

9 * 14 196Ai4 ' 

„ ; 7A 7 8^13^8 

7 1 14 196^14 i 

6 1 14 196Ai4 ' 

D , 10Aj 5A 13 A 5 

1L4 3 4Ai3^ 4l 

3 1 14 196^14 1 

2 1 14 196Ai4 ' 

1 196Ai4 J 
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-{12A12 


14A 14 B n 
B13 


- (13Ai 3 


14A14B 

S13 


12 ^ B12 
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14Ai 4 Bio 

B13 


- (13Aia 


UA 14 B 
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12 -j Bn 
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-{WA l0 
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B13 


- (13Al3 
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Bis j Bi 3 j 
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~{9Ag - 


UA 14 B 8 


(13Al3 - 


14A 14 B 12 


B9 , 
Sis' 


Bis 


Bis ' 
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-{8A S - 


UA14B7 
Bis 
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B s 
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^} 

B13 


B13 


B13 
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